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ABSTRACT 

In the presence of rotation-induced mixing, element diffusion still occurs in 
stellar radiative zones, although at a slower rate than in the case of a complete 
stability of the stellar gas. As a consequence, helium settling leads to vertical 
/^-gradients which, due to the meridional circulation, turn into horizontal fluctu- 
ations. Up to now, the feed-back effect of this process on the rotation-induced 
mixing was currently neglected in the computations of abundance variations in 
stellar surfaces, or artificially reduced. Here we analyse its consequences and 
derive an approximate analytical solution in a quasi-stationary case. We also 
discuss the relative importance of the various physical effects which influence the 
meridional circulation velocity. In a second paper (Theado and Vauclair 2002a), 
we will present a complete 2D numerical simulation of this process while a third 
paper (Theado and Vauclair 2002b) will be devoted to special applications to 
Pop I stars. 

Subject headings: stars: abundances; stars: rot at ion- induced mixing; hydrody- 
namics; diffusion 

1. Introduction 

The fundamental importance of element settling (also called microscopic diffusion) in- 
side stars was recognized by the pioneers of stellar structure computations (see references 
in Vauclair and Vauclair 1982). A star is a self-gravitational gaseous sphere, composed of 
all possible kinds of different chemical elements, with different masses and atomic spectra. 
Due to the pressure and thermal gradients and to the selective radiative transfer, individ- 
ual elements diffuse inside the star, one with respect to the other, leading to a slow but 
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effective restructuration. Mixing processes like convection, turbulence, rotation-induced in- 
ternal motions, compete with microscopic diffusion and decrease its efficiency according to 
the strength of the resulting homogeneization. When Von Zeipel (1924) demonstrated that, 
due to centrifugal effects on the equipotentials of gravity, rotating stars could not be in 
radiative equilibrium, he concluded that the resulting thermal imbalance lead to meridional 
circulation which could be the reason why the effects of microscopic diffusion were not seen 
in most stars. 

Much later, when evidences of "chemically peculiar stars" where obtained, astrophysi- 
cists first discussed them in terms of nuclear or spallation reactions. These explanations 
could not however account for the observations (see, for example, Vauclair and Reeves 1972). 
Meanwhile, the idea that these abundance anomalies could be due to microscopic diffusion 
emerged again and gave interesting results (Michaud 1970, Michaud et al. 1976 (MCV 2 ) , 
Vauclair et al. 1978 a and b (V 2 M and V 2 SM)). From then on, for more than two decades, 
this became "the diffusion hypothesis" which was supposed to be added to the physics of 
peculiar stars, for the aim of explaining their anomalies. This interpretation was misleading 
in the sense that it presented element diffusion as a somewhat "ad hoc" process added in 
some cases, forgetting its fundamental character in the structure of self-gravitational spheres. 

The importance of element settling in the stellar evolution process has recently been 
proved in the Sun by helioseismic investigations. It was known for a long time that helium 
and metals should have diffused by about 20% down from the solar convective zone since the 
birth of the Sun up to know (Aller and Chapman 1960). Recent comparisons between the 
sound velocity inside the Sun as computed in the models and as deduced from the inversion 
of seismic modes, confirmed that this settling really occurred. 

In fact, element diffusion is always present in stars, unless hydrodynamical processes 
increase its time scale in such a way that it becomes much larger than the stellar lifetime. 
In this framework, the basic difficulty is not to account for the occurrence of abundance 
anomalies in some stars, but to understand why the consequences of diffusion processes are 
not seen in all the stars and also why, in chemically peculiar stars, the observed anomalies 
are not as strong as predicted by the theory. This is related to the hydrodynamical processes 
which occur in stellar radiative zones and may slow down the element settling. Rotation- 
induced mixing is certainly one of the most important of these processes. 

The process of meridional circulation and its consequences have been studied by many 
authors, including Eddington (1926), Sweet (1950), Mestel (1953, 1957, 1961), Tassoul and 
Tassoul (1982, 1983, 1989), Zahn (1992, 1993), Maeder and Zahn (1998). The way such 
processes can slow down diffusion have also been studied many times in the literature, 
beginning with Schatzman (1977), V 2 SM and V 2 M. 
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Zahn (1992, 1993) suggested that the transport of angular momentum induced by the 
meridional flow could lead to shear flow instabilities. Due to the density stratification, these 
instabilities should create anisotropic turbulence, more important in the horizontal than in 
the vertical direction. The coupling between the meridional advection and the horizontal 
turbulence would lead to a special kind of mixing for the chemical species, parametrized 
as an effective diffusion coefficient. Such a description has been used, for example, by 
Pinsonneault, Deliyannis and Demarque (1992), Charbonnel, Vauclair and Zahn (1992), 
Richard et al (1996) (RVCD). 

Later on, the transport of chemical species and angular momentum were consistently 
coupled in the computations, with the assumption that no other process interferred with 
rotation-induced mixing (Zahn, Talon and Matias 1997, Talon and Charbonnel 1998). The 
results predicted a non negligible differential rotation inside the solar radiative regions, in 
contradiction with helioseismology. The fact that the Sun rotates like a solid body below the 
convection zone proves that angular momentum has to be tranported by another process, 
presumably magnetic fields (Mestel, Moss and Tayler 1988, Charbonneau and Mac Gregor 
1992, Gough and Maclntyre 1998). 

^From the beginning of these computations of rotation-induced mixing, the importance 
of the feed-back effect due nuclearly-induced /i-gradients was well recognized. Mestel and 
Moss (1986) showed how these /z-gradients could slowly stabilize the circulation and expel 
it from the core towards the external layers. They called this process "creeping paralysis" . 
Indeed, without such a stabilization of the meridional circulation, the stars could not become 
red giants in the correct time scale. 

However, in all these computations the feed-back effect due to the diffusion-induced 
/x-gradients on the meridional circulation was not included or artificially reduced to prevent 
numerical instabilities. Vauclair (1999) showed how, in slowly rotating stars, the resulting 
terms in the computations of the circulation velocity could rapidly become of the same order 
as the other terms. Theado and Vauclair (2001) computed these terms for the case of Pop 
II stars and claimed that the induced diffusion-circulation coupling could be the reason for 
the very small dispersion of the lithium abundances observed in halo stars. 

While the \x terms due to element settling have definitely to be taken into account in the 
computations of rotation-induced mixing, the details of the coupling process which occurs 
when these terms become of the same order of magnitude as the classical terms remain 
difficult to handle. In the present paper, we introduce an analytical approach of the process 
in solar type stars and discuss an approximate solution in a quasi-stationary case. In a 
second paper (Theado and Vauclair 2002a, herafter referred to as paper II), we will present 
the results of a 2D numerical simulation which will help visualize the situation and give a 
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prescription for introducing this process in ID stellar evolution codes. The case of earlier 
spectral types, when the so-called "Gratton-Opik term" becomes larger than unity below 
the convective zone, will be discussed in a third paper, with applications to galactic clusters 
(Theado and Vauclair 2002b, paper III). 

Our conclusion will be that the feed-back effect due to diffusion-induced /^-gradients may 
strongly modify the meridional circulation and meanwhile reduce the efficiency of diffusion 
as suggested in Theado and Vauclair 2001. 



2. Rotation-induced mixing 

In rotating stars, due to the centrifugal forces, radiative equilibrium is not satisfied. 
This has to be compensated by a motion of matter with a vertical velocity derived form the 
equation: 

'ds 



P T ^— + u.Vsj = -V.F + pe n (1) 

where e n is the nuclear energy production, negligible in the outer layers, F the radiative flux, 
s the entropy density, p the local density, T the local temperature, u the matter velocity. 

^From eq.(l) we deduce the vertical component of the circulation velocity which, if 
developped on the second Legendre polynomial, may be written (Zahn 1993, Maeder and 
Zahn 1998, herafter MZ98): 

u r = U r P 2 (cos9) (2) 

Following Zahn 1993 and MZ98, all physical parameters are developped on a level surface 
in the form : 

x = x + xP 2 (cos9) (3) 



Then U r is obtained as : 

P 

U r 



(4) 



pgC P T(V ad - V + V^ 

where P is the local pressure, g the gravity, Cp the specific heat at constant pressure, and : 

p lrW 
p 3 g dr 



This parameter was referred to as 6 in MZ98 : we prefer to call it ( to avoid any 
confusion with the latitude angle, specially when we treat 2D numerical simulations (paper 
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II) ; note that, in MZ98 equation (4.38), T was forgotten in the second term of the bracket, 
as well as in the first term of equation (4.36). 

We have neglected the deviations from perfect gas law, which are not important for 
the purpose of the present paper. The thermodynamical parameters are averaged over level 
surfaces : in the following the bars will be omitted for simplicity. Following MZ98, L is the 
luminosity at radius r and M* — M(l — 2 ^ gp ) where M is the stellar mass at radius r and 
p m the mean density inside the sphere of radius r. 



In most stellar situations, the term 



1-Kgpr, 



is negligible. In the present paper, we discuss 



the case of solar type stars where the "Gratton-Opik term 



» n 2 



2wgp 



(larger than the previous 



term as p replaces p m in the denominator) is always smaller than one below the convective 
zones. Thus, for simplicity, we will replace M* by M in the following equations. The case 
where the "Gratton-Opik term" is not negligible in radiative zones (e.g. in F and A stars), 
will be studied in a forthcoming paper (Theado and Vauclair 2002b, paper III). 

In equation (3), the terms related to the /z-gradients are gathered in while the 
classical terms appear in Eq. When the energy production terms are negligible (as in the 
regions below the outer convective zone, when element settling is important), Eq and are 
given by: 



En 



8 /nv\ 

3 V GMj 



1 - 



Q 2 



Pr, 



3 dr \ dr I r 3 



and : 



E„ 



Pm fr d 
p \ 3 dr 



H, 



dA 

dr 



- (Xa + XT + 1)A 



2H T A 



(6) 



(7) 



Here H T is the temperature scale height; A represents the horizontal p fluctuations = ; Xu 

p 

and xt represent the derivatives: 



dlnx 
dlnp 



XT = 



P,T 



dlnx 
d\nT 



(8) 



^From equation (3), we find that three terms compete in the computation of the merid- 
ional circulation velocity : Eq, E^ and the third one, related to the density fluctuations, 
which we will write : 



E 



c 



M 

L ICp di 



(9) 



A fourth term has to be added in case of large horizontal turbulence, parametrized with 
a horizontal diffusion coefficient Dh- It is due to the fact that such a turbulence modifies 
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the radiative energy transport. In MZ98, this term is added in the expression of Eq (their 
equations 4.36 and 4.37). We prefer to add it separately as : 

M 6 

E h = T - 2 C p TD h C (10) 

which can also be written : 

p r K 

where K is the thermal diffusivity (Zahn 1993). In this way, equation (3) may be written : 

P L 

Ur = p 9 Tc P (V ad - V + V,) M Etot (12) 

with : 

E tot = E n + E^ + E c + E h (13) 

Note that the fourth term E h has a different behavior from the other ones, as it is 
formally proportional to U r when the horizontal diffusion coefficient Dh is approximated as 
in MZ98, that is D h = Cyj • U r where Ch is an unknown coefficient. Rigorously, it should 
be introduced in a different way, on the left side of equation (11). We prefer to leave it 
presently on the right side and will show later (section 3.3) how the prescription given for 
its parametrisation induces that it must be negligible compared to the other terms. 

Then the horizontal component of the meridional velocity is deduced from the equation 
of matter conservation: 

div (w r pr 2 ) = —u g pr (14) 
so that : ^ 

uq = — (pr 2 U r ) sin 9 cos 9 (15) 

2pr dr v ' 

In the classical Eddington Sweet circulation, only Eq is introduced in the computations. 
In this case, the divergence of the upward flux is negative : for mass conservation, matter 
flows horizontaly from the upward flux towards the downwards flux. This situation may be 
dramatically modified when the other terms are introduced, and particularly in the situation 
when E^ becomes of the same order as Eq while the two other terms are smaller. 



3. The importance of /i-gradients 

According to the period in stellar evolution, the four terms included in equation (11) 
may have various relative importance. Here we wish to compare their orders of magnitude in 
solar type stars. Before that, we have to discuss the various ways of computing the horizontal 
p, gradients A which enter the E^ term. 



-7- 



3.1. computations of A 

Due to the meridional circulation, vertical \x— gradients give rise to horizontal ones as 
the upward flow brings up matter with a larger /i while the downward flow brings down 
matter with a smaller fi. The resulting horizontal /z-gradients depend on whether horizontal 
turbulence is present or not. 

In case of a laminar circulation (pure advection), where u r and ug are of the same order, 
the order of magnitude of A is simply given by: 

A~r-Vln/i (16) 

(the case where u r and uq have different orders of magnitude will be came upon and studied 
in section 4) 

According to Zahn (1992) and Chaboyer and Zahn (1992), the shears induced by the 
circulation lead to horizontal turbulence, which reduces the horizontal /i-gradients. In this 
case, introducing a horizontal diffusion coefficient D h , the local variations of \i are solutions 
of: 

i + <r)-i = ~> D ^ (17) 
In the stationary case, the horizontal /z-gradient including turbulence is given by : 

With the assumption that : D h ~ U r ■ r, this equation becomes: 

A L ?hE (19) 

6C h dlnr { } 

where C h should be of order unity (MZ98, equation 4.41 ; note that we have chosen here the 
MZ98 definition of Ch, which is the inverse of its definition in Chaboyer and Zahn 1992) 

At the beginning of the stellar evolution on the main sequence, A is equal to zero. In the 
case of negligible differential rotation, it increases and eventually reaches the critical value 
for which ~ \Eq\. Let us call this value A° crit . It only depends on the physical conditions 
in the star, without any arbitrary parameter. If we neglect the derivatives of A (in paper II, 
all the terms will be included in the complete 2D simulation), we find : 

Kru = En^^T (20) 

where Eq represents the part of the expression of Eq which remains when the differential 
rotation is negligible. 
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We see that A° crit varies like VL 2 . In solar type stars, it is of order 10~ 6 for rotation 
velocities of order 5 km.s -1 (Theado and Vauclair 2001). In case of non negligible differential 
rotation, the value of A crit may become larger than A° crit . 

When there is a non negligible radial differential rotation in the stellar radiative zone, 
C > 0. A special situation occurs if ( becomes equal to A. The equation of state leads to : 



or : 



p = -- + a 

P T n 



T 

T =A ~ C 



(21) 



(22) 



(if the equation of state is different from perfect gas law, factors of order unity have to 
be taken into account in these equations ; we neglect them in the present discussion for 
simplicity) 

If C = , then T/p: is constant along a level surface. If ( — A, the temperature 
fluctuations vanish and the circulation stops unless forced by another process like angular 
momentum transport. In this special case we may write Eq as : 



E, 



n 





r n 2 ' 
i 


Pm 


r d 1 


\GMJ 


2ngp_ 


P 


3d^ \ 



' dA \ 2H T . 2, 

H T — - xt A A + -A 

dr I r 6 



(23) 



Then A crit becomes 



A 



3E% 



crit 



(24) 



which is about 100 times larger than A Q crit below the convective zone of solar type stars. 



3.2. comparison between En, E^, E^ and E h 

Suppose that ( increases rapidly from to A cr u at the beginning of the stellar lifetime 
on the main sequence. Writing A crit ~ aE^ we obtain : 



(25) 



As E^ varies like f2 2 , equation (26) may be written : 

MTC v ^ n 1 dVL 



(26) 
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so that may be larger than only if 



1 dfl 2 L „ 

> «TFF7T^ (27) 



f] 2 di MT C p 

In consequence, the term may be dominant in equation (11) during rotational brak- 
ing, if the braking time scale is smaller than : 

MT C v 

At < a — £ (28) 

For solar type stars, this time scale is of order 10 14 s, i.e. a few million years. 

It is clear that the term may be important only at the beginning of the stellar 
lifetime, and only if the time scale for the angular momentum transport is smaller than a 
few million years. Then, it becomes negligible compared to the two other terms, Eq and 
.This situation may be different in late stages of stellar evolution, as the internal structure 
of the stars then evolve in smaller time scales, leading to a more rapid restructuration, while 
diffusion processes are slowed down because of deep convective zone. However /z-barriers 
induced by nuclear reactions then have to be taken into account, which may change the 
landscape. 

Let us now discuss the E h term in equation (13), as given by equation (10). This term is 
proportional to ( and vanishes in case of negligible differential rotation. Suppose that, on the 
contrary, this term becomes preponderant in equation (12). Replacing Dh by its expression 
Cij • U r , this equation becomes : 

U r = _ m P ~ C h rU r { (29) 

P9 (V ad - V + V„) r 2 

U r being on both sides disappears, and the coefficients in front of Ch( are completely deter- 
mined by the stellar structure, so that the product of these two terms should be fixed, which 
does not seem credible. At the bottom of convective zones in solar type stars, the order of 
magnitude of the product Ch( as derived from equation 30 is typically one. As ( is of order 
10~ 6 to 10~ 4 maximum, it would need a coefficient Ch of order 10 4 to 10 6 , much larger than 
one. In consequence, E h is always negligible in the computations. 

In the following section, we discuss the situation which happens when E^ and Eh are 
negligible and \E^\ becomes of the same order as \Eq\. 
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4. Analytical solution in a quasi-stationary case 

Due to helium settling below the stellar convective zones in solar-type stars, a /x-gradient 
can built so that \E^\ becomes very close to \Eq\ in a time scale which depends on the local 
conditions but may be much shorter than the main-sequence lifetime (Vauclair 1999, Theado 
and Vauclair 2001). In the case when « Eq, a local process emerges where too "big 
terms" (E^ and Eq) nearly cancel each other, thereby creating numerical instabilities in 
ordinary "diffusion and mixing" computations. In Vauclair 1999 and Theado and Vauclair 
2001, we suggested that both the mixing and the settling were then strongly modified. Here 
we try to understand in more details what physically happens in the star when this situation 
occurs. We study it first in an approximate analytical way and later, in a second paper, we 
will present a complete 2D numerical simulation of this process (Theado and Vauclair 2002a, 
paper II). 

If is negligible, and E^ nearly cancels Eq, we may expect, at first sight, that rotation- 
induced mixing freezes out, like in the nuclear core of the star (Mestel and Moss 1986). This 
"creeping paralysis" would first occur in the region just below the convective zone, and slowly 
extend down to deeper layers. However such a scheme forgets that helium settling would go 
on in this case, changing both the vertical and horizontal /z-gradients with time. Thus no 
equilibrium could be reached, and circulation would proceed again. This is the situation we 
analyse below. 

Suppose that the circulation is nearly frozen in a region just below the convective zone. 
This means that on every level surface, the horizontal /z-gradient A is just equal to the critical 
value A crit , for which \E^\ is equal to \Eq\. However, helium diffuses out of the convective 
zone where it is completely homogeneous, that is A = 0. Because of this diffusion, the 
horizontal /z-gradient is always forced to remain below the critical value in a boundary layer 
just below the convective zone (in this analysis, we do not take into account any overshooting 
or tachocline, but it would reinforce this conclusion). Under such conditions, the radiative 
unbalance which is responsible for the meridional circulation and subsequent mixing still 
exists in this boundary layer while it is mostly cancelled out in the frozen region below (but 
it still exists deeper in the star). 

Let us go back to the equation of matter conservation (14). In normal (classical) merid- 
ional circulation, the flux divergence is negative and matter is expelled horizontally from 
the upward flow to the downward flow. Here, on the contrary, as the vertical circulation 
velocity becomes very small below the convective zone except in the boundary layers where 
it increases rapidly upwards, the flux divergence is positive and matter is "aspirated" instead 
of being expelled. We expect that this situation must result in a local loop of matter which 
closes horizontaly from the downward flow to the upward flow (Figure 1). This may, in first 
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approximation, lead to a quasi- stationary stage that we discuss below. 




Fig. 1. — Sketch of the circulation streamlines which may occur when the two terms E^. and 
Eft are of the same order of magnitude. Normal meridional circulation goes on in the internal 
regions where the diffusion-induced /i-gradients are small (note that deeper in the star the 
nuclearly-induced /x-gradients should be taken into account, which are not included here). In 
the "frozen region", the circulation is negligible except just below the convective zone where 
a horizontal current is driven due to element settling. Exact streamlines computed from a 
2D numerical simulation will be given in paper II. 

We assume that A is equal to A cr j t below the convective zone except in a few layers 
where it decreases exponentially from A crit to 0. That is : 



where if is a scale height which may be identified to the diffusion scale height and r c the 
radius at the bottom of the convective zone. 

Let us call Uq the value of U r obtained when = (classical meridional circulation). 
When ^ , we may write : 




(30) 




(31) 



Assuming, in first approximation, that E^ oc A, equation (32) becomes: 



(32) 



- 12 - 



We may now calculate, from equation (14), the magnitude Ug of the horizontal velocity 
which we define by : ug = Ug sin 9 cos 9 



Neglecting the p (r) variation in the boundary region, for a first approximation, we 
obtain : 

Ug ^ -U n 



1 + — 
2H 



[ '—Y^) (33) 



as H < r , this equation becomes : 



r 



t^-^exp^-^j (34) 
Thus the horizontal velocity is much larger that the vertical one for matter conservation: 



Ug > U r (35) 

So matter which settles down below the convective zone is swept up again through this 
boundary loop. A quasi-stationary stage may take place for which the horizontal \i— gradient 
A does not change with time, except for the slow increase of the scale height H. 

In this case, the continuity equation may be written in first approximation : 

^ 9 rd9 ^ T dr 



so that A must be written : 



U r .rd\n.n ^ <91n/x 
Ug dr dr 



The vertical \i— gradient is fixed from this relation : 

d In /I A A crit 



dr ~ 2H 2H 



1 — exp 



H 



(38) 



Although these computations are approximate, they show how the coupling between 
helium settling and rotation-induced mixing can slow down both the mixing and the set- 
tling when the /x-terms are of the same order as the f2-terms in the velocity of meridional 
circulation. 

This does not completely stop the diffusion, but matter flows ajust themselves so that, 
in the layers where this process occurs, the vertical /z-gradients remain constant. We see 
from figure 1 that matter is flowing in opposite directions below the convective zone and 
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deeper in the star. This could create shear flow instabilities in-between. However, as will be 
discussed in paper II, the frozen zone is large enough and is the site of a stabilizing vertical 
/x-gradient, so that we expect this effect to be negligible. As diffusion proceeds slowly, the H 
scale height increases with time : it is not a complete stationary stage. But we expect the 
element settling to be significantly reduced by this process, which may completely change the 
landscape of abundance variations in slowly rotating stars, as already suggested by Vauclair 
(1999) and Theado and Vauclair (2001). 



5. Discussion 



In the present paper, we have discussed the orders of magnitude of the various terms 
which enter the computation of the meridional circulation velocity in solar type stars. We 
have shown in an approximate way what happens when the terms due to the diffusion-induced 
//-gradient become of the same order as the classical terms in the circulation velocity. In this 
case the pattern of the meridional circulation is dramatically modified. 

When the /x-terms are not taken into account (classical circulation), a prescription may 
be used, as given by Zahn (1992), to compute the transport of chemical elements. Due to 
angular momentum transport, a shear flow instability may develop which leads to anisotropic 
turbulence, smaller in the vertical than in the horizontal direction because of the density 
gradient. The horizontal turbulence is parametrised by a turbulent diffusion coefficient 
Dh — ChT-U r while the vertical transport, which is a combination of the meridional advection 
and the horizontal turbulence, is treated with an effective diffusion coefficient : 



When the /x-terms are taken into account, the prescription must be modified. Horizontal 
turbulence has two effects : first it slows down the construction of horizontal /i-gradients, 
thereby delaying the time when the /i-terms become of the same order of magnitude as the 
fi-terms ; second it helps transporting energy horizontaly (the so-called Eh term in equation 
(12)). We have shown however that the Eh term can never be preponderant, because of its 
dependance on U r , which is consistent with the fact that C h should be of order one. In this 
discussed in Vauclair and Theado (2001) for halo stars and will be discussed again 
for other cases (paper II and paper III), the /i-terms and the f2-terms become of the same 
order during the main sequence lifetime of most solar type stars. Then the process described 
here has to be taken into account. 
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We have shown that, when \E^\ becomes very close to \En\, matter flows may ajust 
so that the vertical /i-gradient V In \x remains constant in the "frozen zone" . This "self- 
regulating process" will be confirmed by the 2D numerical simulation (paper II). It can be 
used as a prescription in stellar evolution codes for the computations of abundance variations 
of chemical elements. It assumes however that the time scale for restoring the equilibrium 
/^-gradient is short compared to the time scale of the element settling process which tends 
to destabilize it. In other words, it assumes that Uq is larger than the microscopic diffusion 
velocity. This is the case below the convective zone, but it may not be the case deeper in 
the star where Un decreases quicker than the diffusion velocity. This is why, in Vauclair and 
Theado (2001), we assumed that the "self-regulating " process occurred only in the region 
where Uq was larger than V^/ and that microscopic diffusion occurred unaltered below. 

The question may arise of the orders of magnitude of the other terms in equation (13) 
compared to the difference \Eq\ — l-E^I when this difference is small. In any case, as we have 
seen, the term Eh can never become preponderant. What about the Eq term ? 

E^ is given by equations (26) and (27). We know, from rotational braking in solar 
type Pop I stars, that the external rotation velocity decreases rapidly (in a few million 
years) from about lOOkm.s -1 to a few km.s^ 1 and then go on decreasing very slowly. The 
situation inside the star depends on the process of angular momentum transport. We know, 
from helioseismology, that a transport process more efficient than the one induced by the 
circulation is needed for the angular momentum. As discussed in section 3.3, may be 
important at the beginning of the stellar lifetime on the main-sequence, which is consistent 
with the rotation braking time scale. For it to be also important later, during the self- 
regulating process, it would need a time scale (equation (29) modified) : 



E n MT C„ 

At < " — £ (40) 

En + E^ L 



which can also be written 



Ai<(1 __^,-^ (41) 

^crit Li 

With At of order 10 17 s ( a few billion years), equation (41) can only be satisfied if 



(1 - A) < 10- 3 (42) 
which, in our exponential description, may occur only below a radius r given by: 



r c — r ~ 7H 



(43) 



In summary, mixing induced by the transport of angular momentum may be important 
at the beginning of the main sequence stellar lifetime. Later on, the situation below the 
convective zone is modulated by the competition between the so-called f2-currents, due to 
the classical meridional circulation and the so-called /i-currents, due to the diffusion-induced 
//-gradient. In the present description, we have not taken into account other mixing processes 
which can decrease the horizontal chemical fluctuations, like tachocline layers. This should be 
added in future studies. As proposed in Theado and Vauclair (2001), the reason why no large 
lithium depletion is observed in most of Pop II main-sequence stars, in contradiction with 
Pop I stars, could be due to a different rotation history : as the two stellar populations were 
formed in the Galaxy in different circumstances and at different ages, it could be that Pop II 
stars suffered less rotation braking than Pop I stars so that the final lithium abundance and 
the very small dispersion observed would be a consequence of the "self-regulating process" . 

This process, described here in an approximate physical way, will be studied in paper 
II with a 2D numerical simulation. 

We thank Jean-Paul Zahn for fruitful discussions, critics and comments. 
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